Bose-Einstein condensates on tilted lattices: coherent, chaotic and subdiffusive 

dynamics 
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The dynamics of a (quasi) one-dimensional interacting atomic Bose-Einstein condensate in a tilted 
optical lattice is studied in a discrete mean- field approximation, i.e., in terms of the discrete non- 
linear Schrodinger equation. If the static field is varied the system shows a plethora of dynamical 
phenomena. In the strong field limit we demonstrate the existence of (almost) non-spreading states 
which remain localized on the lattice region populated initially and show coherent Bloch oscillations 
with fractional revivals in the momentum space (so called quantum carpets). With decreasing field, 
the dynamics becomes irregular, however, still confined in configuration space. For even weaker 
fields we find sub-diffusive dynamics with a wave-packet width spreading as t^^^. 

PACS numbers: 03.75.Kk, 03.75.Nt 



Recently much attention has been payed to Bloch os- 
cillations (BO) of a Bose- Einstein condensate (BEG) of 
cold atoms in tilted optical lattices (see Refs. [H, [j, 0, 
m, 0, 0, B S [lO, [Hi, to cite few of the relevant papers). 
The main question one addresses here, both experimen- 
tally and theoretically, is the effect of atom-atom interac- 
tions on the otherwise perfectly periodic atomic dynam- 
ics. Although all experiments on BEC's BO were done for 
a localized initial conditions, theoretically a simpler case 
of uniform initial conditions, where atoms are delocalized 
over whole lattice, is usually analyzed [5, ^, [HI, [ll| . It was 
found that depending on the system parameters a BEG 
of interacting atoms may show three different regimes of 
BO [11]: (i) exponentially decaying BO where the con- 
densate rapidly decoheres, which corresponds to unstable 
(chaotic) BO in the mean-field approach; (ii) persistent 
BO and (iii) modulated BO, which both correspond to 
stable dynamics in the mean-field approach. Assuming 
the magnitude of a static field to be the control param- 
eter, these three regimes refer to weak, moderate, and 
strong static fields, respectively. 

In the present work we extend previous studies to the 
case of nonuniform initial conditions, realized in a labora- 
tory experiment. In what follows we focus on the mean- 
field analysis, - a comparison with quantum mechanical 
treatment remains an open and challenging problem, as 
it will be discussed in some more detail in the concluding 
paragraph of the paper. 

In the mean-field approximation, the dynamics of a 
BEG in a tilted optical lattice can be described by the 
discrete nonlinear Schrodinger equation (DNLSE) 



i hdi = dFlai 
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where ai is the BEG complex amplitude in the Ith po- 
tential well. The static force F leads to a linear increase 
of the onsite energy dFl {d is the lattice period), J mea- 
sures the tunneling transitions between the wells and g 
characterizes the atomic interaction. Equation ([T]) ap- 



pears as a canonical equation of motion i Tidi = dH/dal^ 
ihdi = —dH/dai of a 'classical' Hamiltonian function 
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|a^p is conserved under the evolution and 
|a^p = 1 in the following (which assumes 
that g is proportional the number of atoms), as well as 
units where h = 1 and d = 1. 

The gauge transformation ai{t) ex.-p{—iFlt)ai{t) 
converts the static term in Eq. ([1]) into a periodic driving 
with the Bloch frequency F, 
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and a Fourier transform of the amplitudes ai yields the 
Bloch-waves representation bk = Xl^i exp(— m/) a/, 
where tz = 27rk/L is the quasimomentum (— tt < k < tt). 
As follows from (|3]), the amplitudes obey the equation 
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where S{k) is the Kronecker function modulo L. It is 
easy to see that the displayed equations of motion can be 
solved in closed form in the case of uniform initial condi- 
tions, ai{0) = const, yielding bo{t) ~ exp (z-^ sin(Ft) — 
ij^t) and bk^o{t) = , i.e. ordinary BO. The stability of 
this solution with respect to weak perturbations has been 
analyzed in Ref. [J, o\ , also in comparison with related 
phenomena in the many-particle Bose-Hubbard model for 
a small number of particles and lattice sites [11]. 

Motivated by recent experiments [10], where a BEG of 
Gesium atoms was prepared in a harmonic trapping po- 
tential in the Thomas-Fermi regime, we chose a Thomas- 
Fermi distribution 



(0) = V/^ - 
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for aP < P and zero otherwise. We shall use a = 0.001, 
which populates lattice sites between / = — 9 and / = +9 
(the value of P is fixed by normalization) . The results of a 
numerical solution of the nonlinear coupled equations ([3|) 
with initial conditions ([5]) and parameters J = 1, ^ = 10 
are presented in Figs. [T] and [2] as a function of time mea- 
sured in units of Trev as defined in Eq. (|8]) below. The 
site populations Pi{t) = \ai{t)\'^ are shown in the upper 
panel and the quasimomentum distributions in the 
lower panels. When the force F is varied, one observers a 
characteristic transition from a highly coherent evolution 
for strong fields to an irregular chaotic motion for inter- 
mediate fields and a nonlinear diffusion for weak fields. 
In the following sections we will discuss the dynamics in 
the three different regimes in some detail. 

Coherent evolution. For strong fields, as in the case 
F = 100 shown in Fig. (TJ the lattice population Pi{t) = 
|a/(t)p is almost frozen without observable structure. On 
the contrary, the quasimomentum distribution shows a 
highly organized pattern, also denoted as a quantum car- 
pet. The initial distribution is reconstructed at time Trev, 
the revival time, and at rational ratios of t/Trev we find 
fractional revivals. Such highly organized patterns have 
been observed recently experimentally for BECs in tilted 
optical lattices [10]. 

The origin of these quantum carpets lies in the initial 
distribution ([5]). Indeed, in the limit F ^ oo the site 
populations are constant [TI], \ai{t)\ ^ \ai{0)\, and the 
phases evolve linearly in time: 



ai{t)^ai{0)exp [-ig\ai{0)\^ t] 
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FIG. 1: Lattice site populations Pi{t) — \ai{t)\^ (upper panel) 
and the quasimomentum distribution (lower panel) as a func- 
tion of time in unites of the revival time Trev = 2ix j got for a 
force F — 100 (parameters ^ = 10, J —V). The site popula- 
tions are frozen and the quasimomentum distribution show a 
highly organized quantum carpet structure. 




I 



0.2 



0.4 0.6 

t/T 



0.8 



FIG. 2: Same as Fig. [T] however for a weaker force F — 10. 
The site populations are still localized but show fluctuations. 
The quasimomentum distribution appears to be irregular af- 
ter a transient time. 



(see also 0). In the Bloch-wave representation we find 



bk ~ ex.p{—ig/3t) exp [i{f<il-\-gal'^t) 
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an expression already given in [10] • The most prominent 
feature of such Gaussian sums, which have been analyzed 
in much detail during recent years (see, e.g., [TsL fl3, 
[l5l Il6l] ) are the revivals and fractional revivals. From 
Eq. ([7j) the revival time is given by 



27r /ga 



(8) 



in agreement with the numerical results shown in Fig. [T] 
for F = 100. In comparison with the Bloch period Tb = 
27r/T we have Trev/Te = F/ga = 10^, i.e. the revival 
time is much larger than the Bloch period. For rational 
fractions t/Trev = m/n {n^m integer), there appear n 
(approximate) copies of the initial distribution of reduced 
size. For irrational ratios we find fractal distributions 

It should, however, be pointed out that the limiting 
strong field behavior in Eqs. (|6]) and ([7j) describes the 
evolution for finite values of F only approximately and 
the coherent evolution shown in Fig. [T] is a transient phe- 
nomenon for intermediate forces. For F = 10, as shown 
in Fig. O the quantum carpets gradually disappear after 
a time interval T^oh 0.2 T^ev 

for F = 10. This coher- 
ence time was found to vary approximately linear with 
F. 

Chaotic evolution. The distortion of the quantum car- 
pets in the quasimomentum distribution for t > Tcoh is 
accompanied by fiuctuations in the site populations Pi{t). 
However, lattice populations are still entirely localized on 
the initial interval and fiuctuate around the initial popu- 
lation Pi{0). A measure for these fiuctuations is given by 
the quantity C{t) = |P^+i(t) — P^(t)|, which saturates 
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FIG. 3: Site populations Pi{t) (upper panel) and dispersion 
M(t) (lower panel) as a function of time in units of the Bloch 
period Tb = 27v/F for F = 0.25. The other parameters are 
the same as in Figs. [1] [2] except the upper limit of the time 
axis, which is four times larger. 

in time on a level, which increases strongly with decreas- 
ing F. This behavior of C{t) may serves as an indicator 
the transition from regular to irregular dynamics. 

An alternative indicator of this transition is the finite 
time Lyapunov exponent 

Xit)=\n\Sa{t)\/t, (9) 

where S3.{t) evolves in the tangent space according to the 
linear equation 

i^=M[a{t)]5a{t), (10) 

where Ai is the Jacobian matrix of the evolution equation 
(|3]). We have found that for moderate and weak fields, 
the finite time Lyapunov exponent converges to a posi- 
tive value after approximately ten Bloch periods, which 
indicates a chaotic motion. For strong fields F > 10 the 
finite time Lyapunov exponent behaves as 1 /t during the 
whole simulation time period [17]. 

Nonlinear diffusion. As mentioned above, for a weak 
field the dynamics of the system is chaotic with no sign 
of the quantum carpet in the quasimomentum distribu- 
tion and erratic evolution of the site populations. This 
regime is depicted in the upper panel of Fig.[3l where one 
notices an additional effect not present in Figs. 1,2- the 
wave packet spreading. More quantitatively this can be 
seen in the lower panel, which displays the dispersion 
M{t) = T.il'^Pi - {Y^i^Pif of the population distribu- 
tion. A closer analysis of the functional dependence of 
M(t) reveals a V^-law with the prefactor increasing with 
decrease of the static field magnitude. In addition to 
Fig. [3] and for purposes of future reference. Fig. [H shows 
the distribution of the site populations Pi at the end of 
the numerical simulations where, to reduce fluctuations 
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FIG. 4: Distributions of site populations Pi at the end of 
the numerical simulation discussed in Fig. [3] (blue solid line) 
compared to the initial distribution (red solid line) on linear 
(left panel) and logarithmic (right panel) scales. 
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FIG. 5: Solution of the nonlinear diffusion equation for D — 
50 on linear (left panel) and logarithmic (right panel) scales. 
Distributions at t = (red solid lines), t/27v — 100 (blue 
dashed lines), and t/27r = 1000 (blue solid lines) are shown. 

of Pi{t)^ they are averaged over the last 25 Bloch periods. 

A qualitative explanation for the observed sub- 
diffusive spreading is as follows. In view of the appar- 
ently random behavior of Pi {t) we replace in Eq. (j3|) the 
squared amplitudes \ai{t)\'^ in the interaction term by a 
random variable, |a^(t)p ~ ^(t), with an exponentially 
decaying correlation function 

R{t - t') = mW) = ?exp(-|t - t'\lT) . 



Assuming the quantity ~ \ai{t)\'^ in the last equation 
to be independent of this leads to a diffusive spread- 
ing of the distribution according to the discrete diffusion 
equation f\3\ 

Pi = D{Pi+i-2Pi+Pi-,) (11) 

where the diffusion coefficient D is given by 

D = 7/(^2 + 7^) ^7/^2 with 7 = ^r. (12) 

(Note that Eq. (pT]) is a discretization of the continuous 
diffusion equation dP/dt = Dd^P/dx^ and conserves the 
norm ^iPi = 1.) Since in the present case the quantities 
\ai{t)\^ = Pi{t) depend on / , it is plausible to assume that 
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the diflFusion coefficient locally depends on / as = 
DP^{t). Then the site populations Pi obey a nonlinear 
diffusion equation 

Pi = D (Pf+i - 2Pf + . (13) 

This equation can be considered as a discretization of the 
continuous nonlinear diffusion equation 

dP - d dP 

which appears, for example, in the problem of gas diffu- 
sion in a porous media [19]. For the considered case u = 2 
one of the exact solutions of Eq. (p!4|) is a semicircular 
distribution with a radius growing as t^/^, which implies 
that the second momentum increases as M{t) ~ 

The discrete nonlinear diffusion equation (p!3|) seems to 
inherit properties of its continuous counterpart, although 
we are not aware of any formal proof of this statement. 
As an example. Fig. [5] shows the result of numerical so- 
lution of Eq. (p!3|) for the initial Thomas- Fermi distribu- 
tion. Note, that asymptotically the solution is insensitive 
to the particular shape of the initial distribution and one 
gets the same result, for example, for a Gaussian of the 
same width. 

A comparison between Fig. [H and Fig. [5] shows that 
there are deviations of the actual profile for the site pop- 
ulations from that predicted by the nonlinear diffusion 
model (p!3|) . Nevertheless, the model (p!3|) is capable to 
capture some essential features of the system dynamics, 
in particular, the t^/^-law for the wave-packet spreading. 

In conclusion, we have considered Bloch dynamics of 
a BEG in tilted optical lattices for the Thomas-Fermi 
initial profile, a typical initial condition in a laboratory 
experiment. Similar to the case of uniform initial con- 
ditions studied earlier, the system dynamics is found to 
strongly depend on the magnitude of a static field: it is 



regular in the strong field limit and chaotic in the weak 
field limit. In the former case of regular dynamics the 
BEG wave packet is frozen in the configuration space 
and shows quantum carpet evolution in the momentum 
space, a regime already observed in the experiment |Tq| . 
In the latter case of chaotic dynamics the BEG spreads 
sub-diffusively, with the second momentum growing ap- 
proximately as V^. This latter regime deserves further 
studies because of its relation to the problem of quan- 
tum chaotic diffusion. Indeed, it has been known since 
the pioneering work [20] that quantum interference ef- 
fects strongly modify the classical diffusion due to chaotic 
dynamics (the so-called phenomenon of dynamical local- 
ization [ill]). Because the mean- field approach used in 
this work can be considered as a pseudo classical limit of 
the quantum many-body dynamics [11], it is very inter- 
esting to compare the above mean-field prediction with 
the behavior of the second momentum calculated on the 
basis of the Bose-Hubbard model, as well as that ob- 
served in a laboratory experiment. 

Subsequent to the preparation of this paper we learned 
about recent work [22] where the authors address the 
same problem of different dynamical regimes in biased 
DNLSE, although with a different motivation. As ini- 
tial conditions an extreme case of single site population 
was mostly considered. Among other regimes, a regime 
with no wave-packet spreading and one with sub-diffusive 
spreading, where M{t) ~ t^-^^, were reported. Taking 
into account that solutions of a nonlinear equation are 
typically sensitive (even asymptotically) to the type of 
initial conditions, results of the both papers may be con- 
sidered as consistent. 
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